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ABSTRACT 

Aims. We aim to characterize the properties of the prestellar and protostellar condensations to understand the star 
formation processes at work in a young HII region 

Methods. We have obtained maps of the 1.25mm thermal dust emission and the molecular gas emission over a region of 
20' X 10' around the Trifid Nebula (M20), with the TRAM 30m and the CSO telescopes as well as in the mid- infrared 
wavelength with ISO and SPITZER. Our survey is sensitive to features down to N(H2) ~ 10^^ cm~^ in column density. 
Results. The cloud material is distributed in fragmented dense gas filaments (n(H2) of a few 10'^ cm~'^) with sizes 
ranging from 1 to 10 pc. A massive filament, WF, with properties typical of Infra Red Dark Clouds, connects M20 
to the W28 supernova remnant. We find that these filaments pre-exist the formation of the Trifid and were originally 
self-gravitating. The fragments produced are very massive (typically 100 Mq or more) and are the progenitors of the 
cometary globules observed at the border of the HII region. We could identify 33 cores, 16 of which are currently forming 
stars. Most of the starless cores have typical H2 densities of a few 10* cm~'^. They are usually gravitationally unbound 
and have low masses of a few Mq. The densest starless cores (several 10^ cm~"^) are located in condensation TCO, 
currently hit by the ionization front, and may be the site for the next generation of stars. The physical gas and dust 
properties of the cometary globules have been studied in detail and have been found very similar. They all are forming 
stars. Several intermediate-mass protostars have been detected in the cometary globules and in the deeply embedded 
cores. Evidence of clustering has been found in the shocked massive cores TC3-TC4-TC5. 

Conclusions. M20 is a good example of massive-star forming region in a turbulent, filamentary molecular cloud. 
Photoionization appears to play a minor role in the formation of the cores. The observed fragmentation is well ex- 
plained by MHD-driven instabilities and is usually not related to M20. We propose that the nearby supernova remnant 
W28 could have triggered the formation of protostellar clusters in nearby dense cores of the Trifid. 

Key words. ISM: clouds — dust, extinction — HII regions — ISM: individual objects : Trifid — Stars: formation 



^ 1. Introduction (Bertoldi, 1989; Lefloch & Lazareff 1994, 1995). Systematic 

studies show that these "stellar factories" tend to form 
_ Many observations suggest that HII regions could play an 1^0^ massive and luminous (proto)stars, with luminosities 
Cj important role in spreading star formation throughout the > iqS (Sugitani et al. 1991, Dobashi 2001). It is es- 
Galaxy. As reviewed in Elmegreen (2002), most of the em- timated that more than 15% of the massive stars form in 
bedded young stellar clusters in the solar neighbourhood cGs. The question of whether the mechanical and radiative 
are adjacent to HII regions excited by more evolved stars : luminosities of OB stars can trigger the formation of sub- 
Orion, Perseus 0B2, the Rosette nebula, W3-5, M17, to sequent generations of stars, in particular of massive stars, 
name but a few. Yamaguchi et al. (1999) estimate that 10%- in these condensations is still an open question. 
30% of stars in mass in the Galaxy are formed under the r^^^ u^^Ug^^ collapse" model, first proposed by 
mfluence of adjacent HII regions. In particular, the con- Elmegreen & Lada (1977), is particularly interesting as it 
densations encountered at the border of the HII regions, gj^^^les the formation of massive condensations and sub- 
the so-called bright-rimmed Globules and cometary glob- gequently of massive objects out of an initially uniform 
ules (CGs), have long been identified as stellar nurseries i^edium. Systematic studies have shown that this mech- 
or "incubators" (Sugitam et al. 1991,1994; Lefloch et al. ^^-^^ ^^rk in some isolated Galactic HII regions 
1997; Dobashi et al. 2001). These condensations are lo- ^^^^ ^^^^^^ responsible for the massive YSOs observed 
cal clumps which emerge from the expanding molecular ^-^^ border of the nebulae (Deharveng et al. 2003,2006; 
layer surrounding the nebula. Several theoretical and ob- Zavagno et al. 2007) . In the case of a turbulent environment, 
servational studies have shown the important role played molecular gas is no longer homogeneous but distributed 
by the Ly-contmuum radiation of the nebula m the evolu- ^j^^py filamentary structures; it is not clear however 
tion of these condensations, which is now well understood ^^at are the most efficient processes that account for the 

formation of the subsequent generations of stars. Several 
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second-generation star formation, and proposed two main 
mechanisms of induced star formation, either direct trig- 
gering from HII regions and other pressure mechanisms 

associated with hot stars, or gravity-driven streaming of 
gas along the filaments to build dense cores (Nakamura, 
Hanawa & Nakano, 1993; Tomisaka, 1995; Fiege & Pudritz 
2000; Fukuda & Hanawa, 2000). 

We have started a comprehensive study of a young HII 
region, the Trifid Nebula (M20), small enough to allow a 
complete census of the young stellar and protostellar pop- 
ulation, to address this question. The Trifid Nebula ap- 
pears as a small dusty nebula in a young evolutionary age, 
at an hehocentric distance of 1.68 kpc (Lynds et al. 1985). 
Mapping of the thermal dust and molecular gas emission 
of the nebula at millimeter wavelengths revealed the pres- 
ence of several massive protostellar cores (dubbed TCI to 
TC4) , some of which are exposed to the ionizing radiation of 
the exciting star (Cernicharo et al. 1998, CL98). ISOCAM 
observations confirmed the extreme youth of the protostel- 
lar sources detected in the cores TC3 and TC4 (see also 
Lefloch & Cernicharo 2000, hereafter LCOO). So far, only 
the cometary globule TC2 has been analyzed in detail from 
its line and continuum emission at mid-infrared, millimeter 
and centimeter wavelengths (Lefloch et al. 2002, herafter 
L02); its protostellar content could not be characterized 
from the ISOCAM data, however. 

The last phases of protostellar evolution in the nebula 

have been identified thanks to ISOCAM and SPITZER. 
Deep images obtained by SPITZER in the range 3.5 — 24/xm 
with IRAC and MIPS have revealed the population of low- 
mass objects which formed together with the exciting star 
of the nebula (Rho et al. 2006). In addition to 32 protostel- 
lar candidates, about 120 yoimg stellar objects (Class II) 
have been detected, which are surrounded by warm dust, 
presumably distributed in protoplanetary disks. Analyzing 
the mid-infrared and centimeter emission of the central stel- 
lar cluster members, we found that their disks (components 
B, C, D) imdergo copious photoionization (Lefloch et al. 
2001, hereafter LOI; Yusef-Zadeh et al. 2005). 

In this article, we report on a large-scale survey of the 
cold dust emission around the Trifid Nebula, complemented 
with molecular line emission data. Our systematic map- 
ping has allowed us to determine the evolutionary status 
of the SPITZER Class O/I candidates; it reveals the pres- 
ence of stellar incubators in the dense molecular material 
surrounding the HII region, and the progenitors of the pro- 
tostellar cores. This allows us to draw a full evolutionary 
sequence for the incubators in the Trifid Nebula, from the 
early stage when they are embedded in low-density mate- 
rial of the parental cloud, as illustrated by the massive pre- 
Orion cores TC3-TC4, to the late stage when the parental 
cocoon is fully dissipated, leaving the protoplanetary disks 
exposed to the ionizing radiation of the O star, as observed 
by LOI and Yusef-Zadeh et al (2005). 

The article is organized as follows : after presenting the 
observational data (Sect. 2), we summarize the main re- 
sults of the continuum survey (Sect. 3). We then present 
in Sects. 4-5-6 detailed studies of typical cores which have 
formed in the Western filament(TCO), the Eastern lane 
(TCI) and in the South of the nebula (TC2 and TC9). We 
discuss in Sect. 7 the implications on the star formation 
history of the Trifid. The main conclusions are summarized 
in Sect. 8. 



2. The Observations 

2.1. Continuum observations 

The ISO observations of the nebula are fully described in 
L02. Spitzer observations are presented and discussed in 
Rho et al. (2006). In all maps presented here the coordi- 
nates arc given in arcsec offsets with respect to the po- 
sition of HD 164429A, the O star exciting the nebula : 
a = 18'^02'"23.55", 6 = -23°01'51" (J2000.0). 

Observations of the 1.25mm continuum emission in the 
Trifid were carried out in March 1996 and March 1997 at 
the IRAM 30m-tclcscope (Pico Veleta, Spain) using the 
MPIfR 19-channel bolometer array. The passband of the 
bolometer has an equivalent bandwidth ~ 70 GHz and is 
centered at 240 GHz. The beam size of the telescope is 11". 
The final map was obtained by combining several individ- 
ual fields centered on the; brightest condensations of the 
nebula. Each field was scanned in the horizontal direction 
by moving the telescope at a speed of 4" per sec; subse- 
quent scans are displaced by 4". We used a chopping sec- 
ondary at 2Hz with a throw of 30 to 90" depending on the 
structure of the region to be mapped. The dual-beam maps 
were reduced with the IRAM software for bolometer-array 
data (Broguiere, Neri & Sievers, 1995), which uses the EKH 
restoration algorithm (Emerson, Klein & Haslam, 1979). 
Zero-order baselines were subtracted from the data prior 
to restoration. A few maps suffered from strong skynoise 
and were processed using a noise-reduction algorithm that 
removes a correlated noise component from each channel, 
based on the signal seen by the other channels of the array 
(Broguiere et al. 1995). Since the 1.3mm continmmi map 
was obtained by combining several smaller maps taken in 
the dual-beam mapping mode, it is insensitive to struc- 
tures more extended than the typical size of an individual 
bolometer map in azimuth, ^ 4.5'. The possibility of fiux 
loss from dual-beam filtering is inherent with the obser- 
vational method and cannot be discarded a priori for ex- 
tended structures. This point was investigated in detail by 
Motte et al. (1998), who found a negligible impact on com- 
pact continuum cores. We note a good agreement between 
the determinations of column density based on dust and 
molecular gas, also in warm gas regions where molecular 
gas depletion is known to be negligible (see e.g. Bacmann 
et al. 2002). This makes us confident that the flux losses 
resulting from dual-beam filtering are rather small, except 
perhaps in the diffuse cloud in the TCI 5 region. 

Calibration was checked against Mars and secondary 
calibrators at the IRAM 30m telescope. The weather condi- 
tions were good and rather stable during the two observing 
sessions. The opacity was monitored every hour on aver- 
age and we found typical zenith opacities between 0.1 and 
0.35. Pointing was checked every hour as the source tran- 
sits at low elevation at Pico Veleta and corrections were 
always found to be lower than 3". The relative calibra- 
tion was found to be good to within 10% and the abso- 
lute uncertainty is estimated to be about 15%. The final 
map has a size of approximately 20' x 10' and the rms is 
8mJy/ll" beam over the whole map. Hence, it is sensitive 
enough to detect at the 3a level gas-|-dust column densities 
of ~ 10^^ cm~^ per telescope beam. 

In order to outline the weak extended dust components 
in the nebula, we convolved our map with a Gaussian of 20" 
HPBW. The size of the condensations is estimated taking 
into account the deconvolution from the main beam and as- 
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Fig. 1. (left) Mid-infrared continuum emission at 8fj,m around M20 as observed with SPITZER (Rho et ah 2006). The 
dotted line traces the hmit of the field observed with the 1.3mm bolometer array, (right) Thermal dust emission at 
1300^771 around M20. Contour levels are 5, 10, 60, 120, 240, 480 mJy/U" beam. The location of the ClassO/I sources 
identified by SPITZER is marked by squares. 



suming a Gaussian flux distribution. All the results quoted 
in this paper are based on the non-degraded map with 6" re- 
solving power. The continuum source properties are derived 
adopting a dust temperature T = 22 K, a spectral index 
P = 2 and a dust opacity coefficient K250 = O.lcm^g"^. 
These parameters were obtained from the analysis of the 
dust emission in the cores TC2 and TC4 (LCOO, L02). 

2.2. Millimeter line observations 

A region of 10 arcmin size centered on the Trifid was 
mapped at full sampling in the CO J = 2^1,J=l-^0, 
13C0 J = 1 ^ and C^^O J = 1 hues with the IRAM 
30m telescope (Pico Veleta, Spain) in July 1996, using the 
"On-the-Fly" technique. We used both and autocorrelator 
and IMHz filterbanks as spectrometers. The autocorrela- 
tor provided a kinematic resolution of 0.2kms~^ for all 
the transitions. The kinematic resolution provided by the 
filterbanks ranges from 3.5 km for the J = 1 — > transi- 
tions of the CO isotopomers to 1.3 km s""'^ for the J = 2 — > 1 
transitions. An "intermediate" reference position was taken 
10 arcmin East of the Trifid; the molecular emission at 
this position was later observed in position switching mode, 
adopting an Off position checked to be free of emission. 

A large map of 0.5 deg x 0.5 deg was obtained at the 
CSO in the CO J = 3 ^ 2 transition during various observ- 
ing runs between 1998 and 1999, also using the "on-the-fly" 
technique with a reference position free of emission. The 
map was sampled at 20", corresponding to approximately 
one beamwidth. The receiver was connected to an AOS 
which provided a kinematic resolution of 0.4 km s^^. 

Observations of the Trifid in the lines of CO J = 1 ^ 0, 
CS J = 2 1 and J = 3 ^ 2 and HCO+ J =1^0 were 
obtained with the SEST telescope (La Silla, Chile) in 
September 1996 and October 1997. The observations were 
carried out in position switching mode, using the same ref- 
erence position as for the OTF IRAM data. We observed a 
total area of 20' x 20' with a spacing of 20" over most of 
the region. Two strips of 2' x 15' at the eastern and western 
border of the field, where the molecular emission is much 
less, were mapped with a 40" sampling. An AOS was used 
as a spectrometer, which provided a spectral resolution of 
43 kHz (0.11 km s^i at the frequency of the CO J = 1 ^ 
line) . 

We have adopted the main-beam brightness tempera- 
ture scale to express the brightness of the lines studied 
here, with the exception of the ^^CO and ^^CO lines. As it 
appears in our observations, bright ^^CO and ^^CO emis- 
sion extends over regions of several tens of arcmin. At such 
scales, the main-beam temperature is no longer a good ap- 
proximation to the actual line brightness because the er- 
ror beam contribution is no longer negligible. The antenna 
temperature corrected for rearward losses T^ appears as a 
better estimate and we adopt this scale for the main iso- 
topologue transitions of CO. The validity of this approx- 
imation was discussed previously by L02 (see Sect. 7.1). 



The observed molecular lines, their frequency, the telescope 
beamwidth, the main-beam efficiency Beff and the forward 
efficiency Fes are summarized in Table 1. 

The physical conditions in the molecular gas associated 
with the condensations were derived from the millimeter 
transitions of CS and the CO isotopomers. The kinetic tem- 
perature was estimated from the CO J = 2 — > 1 antenna 
temperature, and found to lie in the range 25 — 30 K. In 
some cores, a simple LTE analysis of the ^'^CO and C^^O 
lines could be performed and showed very consistent re- 
sults with the values derived from the ^^CO data. In all 
the condensations, the ratio of the ^^CO J = 1 — > to the 
C^^O J = 1 ^ intensities was found to be very close to 
the canonical abundance ratio (8), within the uncertainties. 
It implies that the C^^O line is optically thin, and the ^'^CO 
either optically thin or moderately optically thick (t ~ 1). 
The very good consistency of the results from the different 
isotopic CO lines expressed in different units(C^^O inten- 
sities are expressed in units of main-beam brightness tem- 
perature, unlike ^^CO and ^'^CO) provides good support to 
the hypothesis made in our analysis. 

In the analysis of the dense g traced by CS, we 

have adopted a lower value for the kinetic temperature 
(20 K), comparable to the cold dust temperature derived 
from the spectral energy distribution of the cores detected 
in the survey (see above). The density and column den- 
sity of the high-density gas were estimated from a sim- 
ple analysis of the CS line emission in the Large- Velocity 
Gradient approximation. The linewidth was obtained from 
a Gaussian fit to the line profiles of the J = 2 ^ 1 and 
J = 3 —^ 2 transitions. We have assumed a standard abun- 
dance [CS]/[H2] = 10-9. 

The dynamics of the molecular cloud gas will be dis- 
cussed in detail on the basis of our molecular line emission 
survey and presented in a forthcoming paper (Lefloch 2008, 
in prep.). 

3. Results 

3.1. Dust Filaments 

The presence of "stellar incubators" in M20 was searched 
for by mapping the distribution of the cold dust emission at 
1.3mm. We have mapped the emission of the parental cloud 
over an area of about 15 arcmin in R. A. by 20 arcmin in dec- 
lination. The final map is presented in Fig. [T] It offers a very 
contrasted picture of the nebula with the best angular reso- 
lution available in our millimeter wavelengths observations 
(11"). All the condensations discovered in the survey are 
distributed in a complex network of filamentary structures. 
Some of these filaments are detected as dark lanes against 
the bright nebular background in the optical. The surface 
layers of the UV-heated filaments in the nebula are detected 
in the mid-infrared with ISOCAM (CL98) and SPITZER 
(Rho et al. 2006). Four main filaments could be identified 
on the front side of M20 (Fig. [T]) : the lanes which trisect 
the nebula (Southern, Eastern, Western lanes) and a long 
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Table 1. Millimeter Line Observations. 
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Table 2. Physical properties of the filaments. 
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a(b) stands for a x lO*". 



filament on the Western side (12 pc), detected in absorp- 
tion in the optical, to which wc will refer as the "Western 
filament" (WF; sec also Lefloch et al. 2003). Observations 
of the molecular gas emission has revealed several other 
filaments of similar length, 10 pc or more, extending far 
beyond M20, both in the Northern and Southern direc- 
tions, on the rear side of the nebula. The TC4 region was 
previously characterized as a compressed shocked layer, 
filamentary-like, of accumulated material Southwest of the 
nebula (LCOO), dubbed "TC4f". Our survey shows that the 
filament stretches far away from M20, to a region of diffuse 
emission. South of the nebula. 

The physical properties of the filaments were derived 
from the millimeter thermal dust emission. The average gas 
column density N(H) was obtained from the average flux 
integrated in the region defined by the intensity contour 
lOmJy/11" beam and is of the order of a few 10^^ cm~^. 
The average density was derived dividing the gas column 
density by the mean filament diameter, estimated from the 
contour lOmJy/11" beam. To ease comparison with mod- 
els, we also computed the mass, mass per unit length of 
each filament, and determined the sound speed from the 
C-'^^O line profiles. The physical properties of the filaments 
are summarized in Table 2. The filaments are rather dense, 
from a few 10'^ cm~^ to about lO^cm"^ in the central re- 
gions. WF is noticeably larger and more massive than the 
other filaments, but the molecular gas density and the ki- 
netic energy density are very similar in all of them. The 
physical properties of WF, such as its size, its mass per 
unit length, compare very well with thoss of the Orion fila- 
ments (Maddalena et al. 1986) and the infrared dark clouds 
(IRDCs; see Beuther et al. 2007 for a review). 

All the filaments display evidence of fragmentation at 
small and large scales. WF is fragmented in a long chain 
of massive condensations which stretches down to the W28 

supernova remnant (SNR), located at a distance of 1.9 kpc, 
similar to M20 (Velazquez et al. 2002). We have estimated 



the mass of each condensation by integrating on the contour 
of millimeter continuum emission at lOmJy/11" beam; we 
obtain masses of several hundreds Mq for each condensa- 
tion (see Tabic 4; also Sect. 7). Evidence of star formation 
was previously found in WF, thanks to the discovery of TC3 
(LCOO). Several water masers were discovered in TC3, one 
of which is associated with a SPITZER protostellar object 
(Healy, Hester & Claussen, 2004). 

The detection of Class O/I protostellar sources in TCI 
and TC8 testifies that star formation is going on in the 
Eastern lane. The Southern Lane is an exception as it 
does not display any hint of star formation activity al- 
though two condensations were discovered in the filament. 
Such non star-forming (pre-stellar) cores play an impor- 
tant role as they bring information on the physical con- 
ditions of the parental cloud, before the onset of star for- 
mation. Unfortunately, detailed molecular observations are 
available only for a few of them (TC0,TC6,TC10,TC16). 

3.2. Dust Cores in M20 

The continuum sources and their physical properties (size, 
mass, column density) are listed in Table 3. The sources 
can be divided into two categories : 

- condensations (or massive fragments) with a typical size 
of 100" (Ipc) and a mass of several hundred Mq. Their 

extent is defined by the contour at lOmJy/11" beam. 

- cores with a typical size of a few 10" (0.1 pc) and a mass 
of a few ten of Mq . They are subunits of the fragment and 
their extent is defined by the contour at half power. 

All the cores detected in our survey have a typical size 
in the range 10" - 20" (0.1 - 0.2 pc). The gas density at the 
continuum emission peak ranges between 10^ cm^'^ and a 
few 10^ cm~^. The mean core density is higher by a factor 
of 10, typically, in the regions where evidence of shock has 
been detected. 
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Table 3. Condensations discovered at millimeter wavelengths near M20. Star- forming condensations are marked by with 
af. 
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6.5(22) 


- 


TC5t 


18:01:57.65 


-23:12:44.6 


140" X 80" 


177 


134/223 


807 


850/1418 


9.9(23) 


WF 


TC6 


18:02:03.37 


-23:08:40.6 


120 X 50 


- 


- 


43 


126/210 
- 


3.6(22) 


WF 


TC6A 


18:02:02.83 


c\n f\f\ c\ 

-23:09:39.2 


11 X 11 


3 


18/30 


- 


2.9(22) 


WF 


TC6B 


18:02:01.94 


-23:09:22.0 


11 X 11 


2.5 


18/30 


- 


- 


2.5(22) 


WF 


TC6C 


18:02:03.36 


-23:08:39.7 


-\ -\ II ■\ -i II 

11 X 11 


2 


18/30 


- 


- 


2.3(22) 


WF 


TC6D 


18:02:03.01 


-23:08:25.7 


11 X 11 


2 


18/30 


- 




2.0(22) 


WF 


TC7t 


18:02:03.50 


-23:09:51.2 


T r* f\l 1 1 f\f\ll 

160 X 100 


17 


39/65 


63 


102/170 


8.8(22) 


WF 


i UoJ 


1 C.nO.Qfl QQ 

io.Uz.oU.oo 


oq.nn.oo 
-Zo.UU.ZZ.Z 


oo X o4 


lo 


1 A /0/1 
14/ Z4 


40 


a a / 1 no 

OO/ luy 


/ .U(zz ) 




TC9t 


18:02:27.25 


-23:03:20.3 


17" X 13" 






4 


7/12 


1.6(22) 




TCIO 


18:02:22.3 


-23:02:56.0 


24" X 9" 


2.5 


5/8 


25 


30/51 


2.3(22) 




TCll 


18:02:20.4 


-23:01:24.0 


44" X 28" 




6 


23/38 


9.0(21) 


s 


TC12t 


18:02:31.66 


-22:56:46.2 


56" X 32" 






9 


76/126 


9.1(22) 




TC13t 


18:02:40.73 


-23:00:30.6 


55" X 35" 






3.5 


65/108 


1.2(22) 


E 


TC14t 


18:02:47.26 


-23:06:13.0 


unrcs. 










4.7(21) 




TC15t 


18:02:34.17 


-23:06:53.0 


21" X 11" 


1.4 




4.7 


18/30 


9.3(21) 




TC16t 


18:02:34.46 


-23:08:06.4 


36" X 26" 


1.8 




9.4 


24/49 


1.2(22) 




TC17 


18:02:16.85 


-23:00:51.8 


11" X 11" 






7 


12/20 


4.4(22) 


W 



^ The virial mass was computed assuming an isothermal sphere with density p oc r"'^ : Mv = 3Ra^/G and uniform density 
Mv = 5Ra^ /G; the velocity dispersions a is derived from C^**0 line observations, or '^'^CO if C^*0 is not detected. 
^ Peak column density; 

^ WF= Western filament; E= Eastern Lane; W= Western Lane; S= Southern Lane 
^ L02. 5 LCOO 



Of the 33 condensations identified in our survey and 

listed in Tabic 1, 16 sources are associated with SPITZER 
protostellar candidates (Rho et al. 2006). Most of the young 
stellar objects identified by SPITZER are not associated 
with large amount of cold molecular gas or cold dust emis- 
sion in our survey (N(H2) > lO^^cm"^). This is the case 
of TC14, for which wc failed to detect any associated con- 
tinuum or line emission and the case of the young stellar 
objects IRS 6 and IRS 7 discovered by LOl; they exhibit 
the 9.7/tim silicate band in emission, a direct evidence of 
the weak amount of parental material still present. We have 
searched for outflow signatures towards the protostar can- 
didates. In most of the cases, confusion from the ambient 
molecular cloud is such that it is not possible to discrim- 
inate the outflow wing emission in tracers as common as 
CO. It is necessary to observe tracers with a more con- 
strastcd emission, hence less abundant, to detect the signa- 
ture of the outflowing gas. Unfortunately, it becomes more 
difficult to detect this emission in low-density gas region. 



such as the cometary globules located on the Eastern side 
of M20. In a few sources, however, we detected the high- 
velocity wings of protostellar outflows in the CS and SiO 
line proflles (TC3,TC4,TC5). We also detected a protostel- 
lar Herbig-Haro jet associated with the protostellar core 
TC2 (CL98; Rosado et al. 1999;L02). 



4. The Western Filament 

Our extended map shows protostellar condensations TC5, 

TC7 in addition to the source TC3, previously detected 
by CL98 and prcstcUar condensations TCOO, TCO, TC6, 
distributed in a chain of fragments regularly spaced along 
WF. The SPITZER image of the region (Rho et al. 2006) 
shows that small clusters of a few protostars have formed in 
the massive gas fragments TC3 and TC5. In this section, wc 
present a detailed analysis of the condensation TCO. The 
properties of TC5 and TC6 are discussed in Appendix A. 
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Table 4. Physical properties of the cores in the parental cloud of the Trifid. 



Name 


Emission Peak 


Tk 


N(H2) 
(cm ^) 


N( CO) 

(cm ^) 


N(CS) 
(cm 


/TT \ ^ 

n[R2) 


Mcore^ 


Location 






(K) 


(cm ^) 






TCOA 


(—238", +294") 


25 


4.9(22) 


5.5(16) 


1.8(13) 


5.5(5) 


49 


WF 


TCOB 


(— 302 , +292 ) 


25 


2.9(22) 


4.0(16) 


1.0(13) 


4.0(5) 


11 


WF 


TCOC 


(—268", +275") 


25 


2.9(22) 


4.5(16) 


1.1(13) 


4.0(5) 


8 


WF 


TCOD 


(—270 ,+252 ) 


25 


2.9(22) 


3.8(16) 


6.0(12) 


4.0(5) 


9 


WF 


TCI 


(+16 ,+33 ) 


30 


8.4(22) 


3.6(22) 


1.9(13) 


3.0(5) 


28 


EL 


1 U2J 


( 1 '7o'' 1 on'M 
(+72 ,— i^U ) 


on 
30 


8.0(2/j 


A n/ 1 f;\ 
4.U(i6j 


1 O / 1 o\ 

i.8(13j 


3.U(0 j 


0'7 
II 


ooutn 




(— z4d , — zzU j 


on 
20 


1 ^/'oo^ 
i.9(2oj 


n/ 1 c;^ 
D.0(i6j 


o n/ 1 o^ 
2.U(13j 


l.D(Dj 


nn 
90 


Wr 


TP4At 






1 4?9'^"> 






U.U ^^O j 


^fk 
oo 


WP 
vv r 


TCS* 


(-360", -640") 


20 


5.4(23) 




1.0(14) 


2.3(6) 


177 


WF 


TC6 


(-280", -440") 


30 


1.3(22) 


1.5(16) 


1.2(13) 


8.0(4) 


2-3 


WF 


TC8 


(+94", +87") 


30 


3.5(22) 


5.1(16) 


1.5(13) 


4.0(5) 


6 


EL 


TC9 


(+50", -90") 


30 


< 1.0(22) 


1.6(16) 




3.0(4f 


4 


South 


TCIO 


(-20", -66") 


25 


1.2(22) 


2.4(16) 


1.5(13) 


2.5(5) 


25 


SL 


TCll 


(-42", +30") 


25 




1.4(16) 




1.0(4)^ 


6 


SL 


TC13 


(+240", +70") 


25 


6.0(21) 


1.0(16) 


6.0(12) 


1.0(5) 


3.5 


EL 



^ determined from dust 
^ determined from CS 
* Lefloch et al.(2008). 
t L02. t LCOO. 



Fig. 2. Millimeter continuum and line emission observed 
towards TCO (contours), superposed on a 8/im image of 
M20 obtained with SPITZER. In the lower left corner of 
each panel is drawn the beam solid angle (HPFW) of the 
observations. For each molecular line, the emission is inte- 
grated between 6.5 and 13.5 km s~^. 



Fig. 4. Molecular line emission observed towards the 
cores OA-OB-OC-OD of condensation TCO at offset posi- 
tion (-238", +294"), (-302",+292"), (-268",+264") and 
(—270", +252"), respectively. Fluxes are expressed in units 
of antenna temperature. The dashed line marks the velocity 
of the main body gas. 



Fig. 3. Thermal dust emission map of TCO at 1.3mm ob- 
served with the IRAM 30m telescope. First contour and 
contour interval are 30 and 15mJy/ll" beam respectively. 



Condensation TCO is located in the Northwest of the 
Trifid at the border of the HII region (Figs. In the op- 
tical, it is detected as a region of large obscuration against 
the diffuse background emission (see e.g. Rho et al. 2006). 
This shows that the fragment is located on the front side 
of the nebula, which is consistent with the velocity of the 
main body gas vi^r = +10.2 km s^^. Dust emission in the 
mid-infrared shows the photo-dissociation region as a bright 
straight bar, with a sharp border at the interface with the 
HII region (Fig. [2]). The molecular gas tracers confirm the 
presence of such a strong gradient in the emission. 

The cold dust emission at 1.3mm reveals a condensa- 
tion with typical dimensions of 100" x 80" (Fig. [S]). The 
morphology is quite complex and cannot be easily fit by 
an ellipse, unlike the other dust condensations identified 
in our survey. The contour at half-power delineates four 
cores, TCOA-OD (Fig. ^. The "brightest" one, TCOA, is 
noticeably larger than the other condensations, by a factor 
of 2 typically. Although the other cores have a comparable 
brightness, their smaller sizes result in much smaller masses 
(typically a factor 4). TCOA-OC-OD are located right ahead 
of the ionization front and the photon-dominated region 
(PDR), as traced by the 8/im PAH band. By contrast, the 



core TCOB is further away from the PDR and appears to 
be better protected from the high-energy radiation. 

The physical conditions in the pre-stellar cores TCOA- 
OB-OC-OD are very similar, with column densities N(CS) ^ 
1.0 X 10^^ cm^^ and gas densities n(H2) ~ 4 — 5 x 10^ cm~'^ 
No evidence for outflow emission was found either in low- 
and high-density gas tracers, like in TC3-TC4, which is con- 
sistent with the absence of protostellar activity (see Fig. [J) . 
Line profiles towards OB, OC, and OD are doubly peaked. In 
addition to the main body gas emission, two velocity com- 
ponents are detected, between +7kms~"'^ and +9.2kms~^, 
and between +llkms^^ and +12.5kms^^. In particular 
the kinematic component at " blueshifted" velocities is de- 
tected in the CS line profiles, either as well separated com- 
ponent or as shoulder to the main body line profile. As 
can be seen in Fig. O both secondary components follow 
closely the border of the PDR excited by M20, and most 
likely trace the signature of the shock driven by the ioniza- 
tion front at the surface of the fragment. The widespread 
CS emission detected all over TCO shows evidence for high 
density gas (n(H2) ~ 10^ cm^'^. 



5. The Eastern Lane 

Three main condensations were identified in the Eastern 
lane : TCI, located close to the central stellar cluster 
(Fig. ED, TC8 and TC13, which both lay further East in 
the lane (Fig. [7]). A fourth, starless, condensation was dis- 
covered in the filament, peaking at (—50", —10"). We do not 
report on its properties in detail as it remains undetected 
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Table 5. Line parameters of the various transitions observed (intensity in K, velocity peak and linewidth in kms~^) 
at the brightness peak of the condensations. The intensity of the C^^O and CS lines is expressed in units of main-beam 
brightness temperature. The symbol "-" implies that the relevant transition has not been observed. NG means that the 
line profile is too complex to be fit by simple Gauss functions. 



Name 


Emission Peak 


CO 
J= 1 ^0 


CO 
J = 2 ^ 1 


CO 
J = 3^2 


c^«o 

J = 1 ^ 


"CO 
J = 1 ^ 


CS 
J = 2 ^ 1 


CS 
J = 3^2 


TCOA 


(-238", +294") 


22.0 


23.7 


12.5 


1.7 


12.4 


2.5 


3.1 




10.5 


10.2 


10.0 


10.4 


10.3 


10.5 


10.5 






3.0 


3.0 


3.0 


1.8 


1.8 


1.8 


1.5 


TCOB 


(-302", +292") 


22.5 


22.5 


16.4 


1.7 


14.7 


2.30 


2.67 




9.9 


10.0 


10.2 


10.6 


10.3 


10.60 


10.55 






2.8 


2.8 


3.4 


0.8 


1.5 


1.0 


0.95 


TCOC 


(— 268 , +275 ) 


20 


22.0 


13.8 


1.5 


11.2 


2.5 


2.75 




10.9 


10.9 


11 


10.6 


10.5 


10.5 


10.5 






1.7 


1.7 


2.0 


1.0 


1.3 


1.2 


1.1 


iOOD 


(—270 , +252 ) 


20.0 


19.0 


15.7 


1.1 


6.8 


1.2 


1.4 






9.8 


9.8 


10.1 


10.2 


9.8 


10.4 


10.6 






3.8 


4.0 


3.5 


1.2 


1.8 


1.4 


0.9 




^^-f-io , -f-oo 


1 Q R 


9^1 A 


1 "i Q 
10. y 


1 zt 


1 9 


o.u 


o.U 






3.2 


3.2 


3.55 


3.7 


3.8 


3.55 


3.85 






2.6 


3.2 


3.2 


1.3 


1.7 


1.0 


1.1 


TP91 


^4-79" I9n"^ 
(T I ^ , — IZU ) 


ol.U 


9Q 7 


18 7 
lo. 1 




lO.U 


o.U 


A A 




7.7 


7.5 


8.2 


7.7 


7.65 


7.7 


7.7 






2.2 


2.4 


2.9 


1.0 


1.6 


1.3 


1.3 


J- V-^O-TV 


^— 94fi" — 99n""l 


NG 


NG 


NG 




0.0 


i-.O 


9 70 










21.6 


21.9 


22.1 


20.7 












1.9 


2.9 


5.8* 


1.9 




( — ioU , — z4d ) 








A Q 


( .4 


1 
z.i 


Z.O 










21.3 


22.1 


22.6 


22.5 












1.2 


2.3 


2.7 


2.6 




(-360 , — o40 ) 




JNG 


JNG 




~ 


0.34 


0.37 














21.1 


21.4 
















3.0 


4.6 


TC6 


(-280", -440") 


NG 


NG 


NG 




3.1 2.1 


0.74 0.65"^ 


0.49 0.35* 














18.9 16.2 


20.0 16.5 


20.0 16.5 














2.8 2.5 


3.0 1.6 


2.8 1.6 


TC8 


(+94", +87") 


21.4 


21.4 


15.7 


1 ..' 


9.6 


1.9 


2.2 




3.8 


3.6 


3.85 


4.1 


3.7 


3.5 


3.8 






2.7 


3.0 


3.0 


1.1 


1.7 


1.7 


1.7 


TC9 


(+50 , -90 ) 


NG 


NG 


NG 


0.40 


3.5 


- 


- 










16.5 


16.5 
















1.1 


1.2 






TCIO 


(-20", -66") 


22.8 


18.9 


16.5 


2.0 


10.0 


2.5 






0.4 


0.4 


0.8 


0.3 


0.3 


0.4 








2.4 


2.8 


2.8 


0.8 


1.2 


1.4 




TCll 


(-42", +30") 


16.1 


25.2 


16.2 


1.0 


2.9 


0.62 






1.1 


1.1 


1.7 


1.2 


1.2 


1.1 








2.0 


2.1 


2.6 


0.7 


1.4 


1.4 




TC13 


(+240", +70") 


20.0 


23.0 


16.5 


< 0.2 


4.2 


1.3 








4.2 


4.5 


4.7 




3.8 


3.8 








2.3 


2.9 


3.0 




1.2 


1.2 





^ L02. ■^LOO. Lefloch ct al. 2008 (in prep) f corrected for beam dilution. 

measured on SiO as CS line profiles are self-absorbed. 
* low resolution spectrum (1 kms~^); contamination from outflow wings and a second gas component on the line of sight. 
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Fig. 5. Emission (contours) of the CO J = 2^1 Unc integrated in the velocity intervals : [7.0; 8.8 km s^^] (left), 
[9.1; 11.1 km s^^] (middle), [11.5; 12.7 km s^^] (right). The emission is superposed on the emission of the main body 
gas (greyscale), integrated in the interval [9.1; 11.1 km s~^]. Intensity contours are respectively : 8, 11, 14,.. 29 K km s~^ 
(left); 20,28,36,40,44Kkms~i (middle), 4,5,7,9,11 K km s^^ (right). 



Fig. 6. Cometary globule TCI as seen in the mid-IR with 
ISOCAM {8-15nm) and SPITZER {3.6fim, A.dfim, b.Sfim, 
8.5/iTO). In the bottom right panel the 1.3mm thermal dust 
emission (black contours) superposed on a greyscale im- 
age of the 13.3 — 16.2/im, as observed with ISOCAM using 
the Circular Variable Filter. Contours of cold dust emis- 
sion are 10,20,40,60 to 160 mJy/11" beam by step of 20 
mJy/11" beam. Contours of mid-IR emission range from 
0.25 to 0.65 mJy/pixel by step of 0.1 mJy/pixel. The arrow 
points to the position of the protostellar source. 



Fig. 7. Millimeter line emission observed in the Eastern 
Lane towards TCI and TC8. The emission is integrated 
between and -|-7kms~^. In the lower left corner of each 
panel is drawn the beam solid angle (HPFW) of the ob- 
servations. The line emission integrated between and 
7 km s~^ is shown in the other panels. For the ^'^CO and the 
^^CO transitions the first contour and contour interval are 
0.1 times the flux peak value. For the CS J = 2 ^ 1 line, 
first contour and contour interval are 0.5 K km s^^. For the 
C^^O emission, the first contour and contour interval are 
0.2 times the flux peak value. The positions observed at 
IRAM are indicated (open squares) in the upper left panel. 



Fig. 8. Spatial distribution of the UIBs at 6.2, 7.7, 11.3^m, 
and the [Nell] line in the TC 1 region (contours) superposed 
on an optical Ha image of the region (HST). First contour 
and contour interval are 0.04 Jy/px and 0.03 Jy/px respec- 
tively. 



in our millimeter continuum survey; we mention it for com- 
pleteness. TC8 and TC13 are discussed in Appendix B. We 
discuss in greater detail the properties of TCI, which is the 
only cometary globule of the sample for which mid-infrared 
spectroscopy could be obtained. 

TC 1 appears as a bright condensation in the Eastern 
lane, 40" East of HD 164492 A (see Figs. MH, both in 
the millimeter continuum and in the molecular gas emis- 
sion. The condensation is well separated from the dif- 
fuse dust lane emission. The contour at half power de- 
fines a central dense core of diameter of 0.16 pc and typ- 
ical density n(H2) = 1.8 x lO^cm"'^. The outer region, 
with a diameter of 0.4 pc is characterized by lower den- 
sities (n(Il2) = 2.5 X 10"* cm^"^). The physical parameters 
of TCI (gas column density, mass, size, density) are very 
similar to those derived for the protostellar condensation 
TC2 in the South of the Trifid (see Table 4 and L02). 

The emission of the optical bright rim on the Southern 
border of TC 1 in the optical image of the Trifid (Fig. [5]), 
as well as the detection of the [Nell] line over the globule 
(Fig. [5]) testifies to the condensation being photoionized. 



The dark appearance of the lanes implies that they lie on 
the front side of the ionized nebula, and hence they are il- 
luminated on the rear side. The CVF imaging of the UIBs 
traces the outer parts of the condensation, which looks like 
an almost continuous ring (Fig. [5]). The low emissivity of the 
central regions can be attributed to absorption by the cold 
core detected at millimeter wavelengths where large dust 
and gas column densities are measured, corresponding to 
Av ~ 145, and Ti2.7^m = 2, adopting the reddening law de- 
rived by Lynds et al. (1985). We have searched for a mid-IR 
signature of the protostar discovered in TCI by integrating 
the ISOCAM flux between 13.3 and 16.2^m, where the ra- 
diated flux is the highest in the CVF band, and the extinc- 
tion, longwards of the silicate absorption range, is the weak- 
est. The images of the globule in the IRAC and ISO LWIO 
filter are displayed in Fig. [SI A point source is detected 
very close (within 4") to the cold dust peak (see arrow in 
Fig. [11). The high-sensitivity SPITZER/IRAC images have 
detected emission from the source at shorter wavelengths. 
It is unresolved in these observations (l"-2"). A SED was 
built using the IRAC colors and the MIPS 24/im data (see 
Rho et al. 2006), which yields an estimate of 570 Lq for the 
source luminosity. This value is typical of an intermediate- 
mass protostar. 

A spectrum of the source between 5 and 16/iTO was ob- 
tained with ISOCAM/CVF after subtracting the emission 
of the condensation averaged over the closest neighbour- 
ing pixels. Instrumental problems have affected the long 
wavelength border of the spectrum, resulting in a spurious 
decrease of the flux beyond 16/im. The spectrum in the 
range 5.0 — 15.5/im is shown in Fig. [5] The spectrum ex- 
hibits a deep absorption due to the silicate band at 9.7/im 
and a rising slope longward of 10/zm, characteristic of em- 
bedded sources. Longward of 15/iTO, the spectrum displays 
another dip, due to absorption by the CO2 ice band, cen- 
tered at 15.2/^771. The 7.7fj,m emission band detected in the 
TCI spectrum is probably due to PAHs excited on the rear 
side of the condensation, as was reported in TC2 (L02). The 
spectrum could be approximately fit by two black-bodies at 
temperatures of 250 K and 140 K modified by dust opacity 
law T^, cx ly^'^, a spectral index value typical in the mid-IR 
regime, surrounded by an absorbing cold layer of column 
density N(II) = 8.0 x 10^^ cm'^, and a total CO2 ice col- 
umn density N(C02) = 3.3 x 10^*cm~^, using integrated 
absorption cross sections measured in the laboratory (see 
Cernicharo et al. 2000). 

In the absence of observations resolving the emission re- 
gion, we have assumed the following sizes : 0.2" (350 AU) 
for the 250 K component, comparable to the diameter of the 
disks detected in the central stellar cluster (see LOl), 0.5" 
(800 AU) for the lower temperature component (140 K). 
Under such hypothesis, we estimate gas column densities 
N(H) = 4.0 X 10^0 cm-2 and N(H) = 3.7 x lO^i cm'^, re- 
spectively. The spectrum is thus consistent with the picture 
of a protostar forming at the center of the TCI core. 
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Fig. 9. Spectrum in the range 5 — 17/im of the protostar 
embedded in TC 1 (thick), after subtracting the ambient 
gas component. The spectrum is fit by two black-bodies 
modified by a dust opacity law cx i^^-^a.t 250 K (dashed) 
and 140 K (dotted) respectively. The total fit is displayed 
in thick. 



Fig. 10. Spectra obtained towards TCI and TC8 in 
the transitions C^^O J = 1 ^ 0, ^^CO J = 1 ^ 0, ^^CO 
J=1^0, J = 2— i-1 and J = 3 — > 2 (from bottom to 
top). Flux is expressed in units of antenna temperature. 
The dashed line marks the velocity of the main body gas. 



Fig. 11. Dense gas emission observed in the lines of CS 
J = 2 ^ 1 and HCO+ J = 1 — > at the emission peak of 
TCI and TC8, in the Eastern dust lane. The SEST spectra 
of TC8 are corrected for the main beam dilution effect. The 
IRAM spectra of TCI are expressed in units of main-beam 
brightness temperature. 



Some residual emission arises in the UIBs at 6.2 and 
7.7/iTO. This is reminiscent of the situation in the core of 
star forming globule TC 2, in which no emission at all is 
detected between 5 and 17/im apart from a weak emission 
peak at the position of the 7.7/im UIB. A simple modelling 
showed that this emission arises from the PDR on the rear 
side of the globule, where photoionization takes place. It 
is strongly absorbed by the material of the globule, except 
near 6.2 and 7.7/im where local minima in the dust absorp- 
tion curve allow some mid-IR photons to escape. Hence the 
ring configuration arises naturally if, just like TC2, the rear 
side of TCI undergoes heavy photoionization, once the sur- 
rounding lower-density material has been photoevaporated. 
From the absorption produced by the cold dust in the PAH 
6.2/iTO band and the [Nell] line, we derive Ti2.7pm ~ 1.7 and 
T6fj,m > 2.5. Those values are consistent with the extinction 
derived from the millimeter continuum observations and an 
opacity law Ti, cx i/^'^. 

A secondary component, detected at visr ~ +1.8 km s^^ 
could trace the kinematic signature of the shock front 
driven by the photoionization of the surface layers. 



6. The Southern region 

Cometary globule TC2 is the first (proto) stellar incubator 
identified in the Trifid thanks to the splendid photoion- 
ized jet HH 399 which propagates out of the head of the 
globule (see Fig. [121 also Cernicharo et al. 1998). The dust 
and molecular line emission of cometary globule TC2, as 
observed in the mid-infrared with ISO and at millimeter 
wavelengths, was discussed in detail by L02. The signature 
of the radiatively-driven implosion of the globule was unam- 
biguously detected in the high-density gas tracers emission 
(CS, HCO"'"); numerical modelling of the globule yielded a 
duration of 0.3 Myr for the photoionization. It comes out 
that the lifetime of the globule exposed to photoionization 
is about 1.9 Myr, long enough to allow the protostar to 
complete the accretion phase. 



Fig. 12. The TC2 region seen in the mid-IR with ISOCAM 
(8— 15/im) and SPITZER (3.6/iTO, 4.5/im, 5.8/im, S.O/im). 
In the bottom right panel the 1.3mm thermal dust emission 
(black contours) superposed on a greyscale image of the ion- 
ized gas emission as seen in the optical with HST. Contours 
of cold dust emission are 10, 20, 40,60 ...160 mJy/U". The 
arrow points to the position of the protostellar source. 



Fig. 13. Spectrum in the range 5— 17//m of TC9, after sub- 
tracting the emission of a reference position. The spectrum 
is fitted by two modified blackbodies with a dust opacity 
cx i^^-^, at temperatures of 5000 K (dashed) and 320 K 
(dotted) respectively. The sum of the two components is 
displayed in a thick line. 



The ISO data were not sensitive enough to detect emis- 
sion from the protostar. Based on the SED obtained with 
ISO in the range 50 — 200 fj,m, L02 estimated a source lumi- 
nosity L 500 Lq , from which they concluded the proto- 
star was an intermediate-mass candidate. Rho et al (2006) 
obtained a similar estimate (590 L0) using an 3.6 — 25/im 
SED constructed from the IRAC and MIPS colors. The 
protostar powering the jet was only tentatively identified, 
based on VLA observations of the free-free continuum. The 
much higher sensitivity of the IRAC images has allowed to 
detect a few sources with colors typical of young stellar ob- 
jects, in the region around cometary globule TC2 (see Rho 
et al. 2006). SPITZER detected only one source with proto- 
stellar colors inside the globule. Interestingly, this source is 
aligned in the HH 399 jet direction, and lies close to the cold 
dust emission peak (Fig. [T^ . It coincides with one of the 
sources marginally detected at the VLA (L02). The source 
detected with SPITZER/IRAC is most likely the protostar 
powering HH 399. 

ISOCAM detected emission between 8 and 15/im from 
a small condensation (13" x 17") in the Northwest of TC2 
at offset position (Fig. [T^ . The condensation is not cen- 
tered on the flux peak, which indicates the presence of 
a second, fainter, source. The SPITZER/IRAC images at 
shorter wavelength indeed unveil two sources located at 
(50.7", -89.2") and (49.6", -97.4"), respectively, lying in 
a common envelope of 8" x 17". The two sources have 
very different IRAC colors. The Northern source is iden- 
tified as a Class I, with IRAC fluxes similar to those mea- 
sured towards TC4A, hence compatible with the birth of 
an intermediate-mass object. The Southern source has typ- 
ical stellar colors (Rho et al. 2006). It could be a bright 
field star, located behind the small dust core, or a young 
star which has emerged from the photoevaporated layers, 
like the EGGs of the Elephant trunks in M16 (Hester et al. 
1996). 

Mid-infrared spectroscopy in the range 5 — 17/im was 
obtained with the CVF. The pixel size was 6". We present 
here a spectrum of the emission averaged over the conden- 
sation (Fig. [T5)l . after subtracting the emission of a refer- 
ence position, chosen 1 arcmin North of TC9. The spectrum 
is dominated by continuum emission decreasing between 5 
and 16.5/im, and a strong absorption in the 9.7/im silicate 
band. This shows evidence for a high temperature region in- 
side the object. The overall emission can be roughly repro- 
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Fig. 14. Mid-infrared (ISOCAM) and millimeter line emis- 
sion towards TC9. The molecular line flux is integrated 
between f5.25 and 17.25 km s~^. Solid circles indicate the 
HPBW beam size, (upper left) : The mid- infrared 8— 
is shown in black contours. Levels range from 160 to 180 by 
10, 200 to 240 by 20, 240 to 360 by 40 MJy/sr. In thin con- 
tours is drawn the 1.3mm emission; contours range from 10 
to 30 by 10 and from 40 to 150 by 15 mJy/ll"beam. We in- 
dicate for each panel, the observed transition, the first con- 
tour and the contour interval, respectively, (upper right) : 
0.2 and O.IK km s'^ (middle left) : 0.4 and 0.2Kkms-i. 
(middle right) : 1.6 and 0.2 K km s""'^. (bottom left) : 15 and 
IKkms"^. (bottom right) : 14.5 and 0.5Kkms^^. 



Fig. 15. Molecular line emission observed towards TC9. 
Fluxes are expressed in units of antenna temperature. The 
dashed line draws the velocity of the main body gas of TC9. 



duced by a black-body law at a temperature of 250 — 400 K 
modified by a dust opacity law t^, cx u^-^ , absorbed by a 
cold dust layer with an extinction Av = 30, which corre- 
sponds to a gas column density N(H2) = 1.6 x 10^^ cm^^, 
adopting the reddening law derived by Lynds et al. (1985) 
in the Trifid (Fig. [T3|) . The uncertainties in such a fit are 
rather large due to the narrow wavelength range observed. 
We found that introducing an additional component trac- 
ing photospheric emission at about 5000 K provides a bet- 
ter fit to the spectrum. Assuming sizes of 1 AU for the 
photospheric emission region and of 100 AU for the warm 
dust component, we estimate a typical gas column density 
of 3.2 X 10^" cm~^ in the warm dust region. The assumed 
sizes are arbitrary, as the only constraint we have is that 
they are less than 1", the size of the IRAC PSF at 3.6^m. 
The ISOCAM spectrum is typical of young stellar objects in 
the Class II phase and consistent with the emission arising 
from a young stellar object surrounded by a protoplanetary 
disk, embedded in a small gas condensation. The evolution- 
ary age of Class II objects is typically several 10^ yr up to 
a few Myr, i.e. it compares well with the age of the HII 
region. It suggests that such sources could have formed in 
the same burst of star formation as the exciting star of 
M20. From the gas column density of absorbing material, 
we estimate a mass of 4 — 8 Mq for the TC9 core. 

Searching for molecular gas emission from TC9, we have 
detected a clump of gas in the millimeter lines of CO and its 
isotopologues (see Fig. [T4|). The lines peak at 16.5 km s~^ 
(Fig. [T5)l . which implies that TC9 is located on the rear side 
of the nebula. The clump is hardly resolved at 1.3mm in the 
^^CO J = 2 ^ 1 transition. The kinetic temperature of the 
molecular gas is difficult to estimate as the condensation 
is unresolved by the telescope beam. The ^^CO J = 1 ^ 
line consists of a narrow line of 7 K superposed on a broad 
pedestal of 4 K. Taking into account the beam dilution 
effect on the core, we find that the ^^CO J = 1 ^ line 
brightness has to be less than 35 K. 



7. Discussion 

7.1. Core Properties 

In many cores, self-gravity is comparable to the kinetic en- 
ergy available, hence plays an important dynamical role in 
the evolution of the cores. This is the case of cometary glob- 
ules TCI, TC2, TC8, embedded cores like TC3A, TC4A, 
TC5, all of them are star-forming condensations. However, 
the non star-forming condensations TCOA, TCOO are also 
found in this class. This suggests they might be on the verge 
of gravitational collapse. 

In most of the pre-stellar condensations (e.g. 
TC0B,0C,4C,6A-D,7), self-gravity is much lower than 
the available kinetic energy. This also applies to the 
evolved cometary globules TC12,13,14,15,16 which have 
already formed protostars, detected with SPITZER. In 
such objects, the dense gas emission is restricted to the 
region right behind the bright rim and we do not find 
evidence of a dense protostellar envelope which would 
be associated with the SPITZER sources, like in e.g. the 
young ClassO/I TC2-3. Most likely, the bulk of accretion 
is over in these cometary globules, as only low column 
densities of material are still available there. 

Overall, the star-forming cores detected in our survey 
have masses of 10 M0 or more (up to 100 — 200 Mq; see 
Table 1). The two cometary globules which have been stud- 
ied in detail, TCI and TC2, exhibit very similar physical 
properties, in terms of core mass, density, temperature and 
turbulence. With an average core density of 1.8 x 10^ cm"'^ 
and sound speed a = 0.42 km s~^, the Jeans length is Aj is 
K, 0.17pc in TCI and TC2, hence approximately the core 
diameter. It is therefore no wonder that only one proto- 
star was detected in the condensations. In other cases, the 
Jeans length Aj is less than one core radius and the Jeans 
mass is a small fraction of the total mass available in the 
condensation. Hence, the cores appear to be dense and mas- 
sive enough to give birth to a few stellar objects. This is 
in agreement with the IRAC maps of TC3, TC4, and TC5 
which unveil clusters of protostars. 

On the contrary, most of the non star-forming (pre- 
stellar) condensations have masses of a few Mq , and self- 
gravity is not large enough to counterbalance the kinetic 
energy of the gas. The only exceptions are TCOA and TCOO 
(49 Mq). We conclude that the cores in TCO and TCOO are 
privileged sites for the formation of the next generation 
of stars. We speculate that the gravitational collapse will 
start as soon as the cores will be compressed by the shock 
preceding the ionization front. The prestellar cores in TCO 
and exhibit physical properties (mass, temperature, den- 
sity) similar to those of the star- forming cores TCI and 
TC2. Hence, we expect that most of the new forming stars 
will be low- (and/or intermediate) mass stars and will form 
inside these fragments. 

7.2. Evolutionary Sequence for the cores 

The cold dust cores detected in our survey are submitted 
to a large variety of physical conditions depending on the 
photoionization in the parent cloud : very deeply embedded 
cores in massive molecular gas fragments (e.g. TC3-4-5); 
cometary globules, i.e. cores surrounded by a dense enve- 
lope exposed to the ionizing nebular radiation (e.g. TCl-2); 
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cores in an advanced stage of photoionization, which led 
them to get rid of their envelope (e.g. TC12-13). 

Undoubtedly, such sources are in a very advanced stage 
of protostcUar evolution and photoionization. Some other 
sources, such as TC9 and TC13, are associated with molec- 
ular material detected in line emission, but remain below 
the detection limit of our millimeter continuum survey. 

Our survey allows us to draw the full evolutionary se- 
quence of cometary globules exposed to the ionizing nebu- 
lar radiation, from their initial stage of formation, as illus- 
trated by the pre-stellar cores embedded in the fragments 
of dense molecular gas TCO and TCOO, to the late stages 
when the bulk of the envelope has been photoevaporated, 
like in TC9. Condensations TCl-2-8, directly exposed to 
the harsch nebular ionizing radiation, are still emerging 
from the low-density layers of the parental molecular cloud. 
This differs from the situation in M16, another well-studied 
HII region, where the erosion of molecular gas condensa- 
tions by the ionizing radiation has shaped "pillows" in the 
parent cloud. 

The free- free emission was observed towards TCI at the 
VLA by Yuscf-Zadeh et al. (2005), who derived an electron 
density of « 2000 cm~^ in the bright rim, i.e. similar to 
the value measured in the bright rim of TC2. Taking into 
account that the other properties are similar, the mass-loss 
rate induced by photoevaporation should be close to the 
value derived at the surface of TC2 : 34Af0Myr~^ The 
lifetime of TCI is therefore « 2Myr, similar to the ex- 
pected lifetime of an O 7.5 III star (Martins et al. 2005). 
The photoevaporation timescale of the condensations TCI 
and TC2 is long enough to allow the central protostellar 
objec;ts to accrcite the bulk of material from their enve- 
lope before being exposed to the ionizing radiation. Hence, 
the incubators located in the densest part of the molecular 
cloud, in the dust lanes of the Trifid, will eventually allow 
their young stellar objects to reach "maturity" (the proto- 
planetary phase) before being exposed to the harsh ionizing 
radiation. 

The fate of TC9, which is located in a less dense re- 
gion on the rear face of the nebula, is somewhat different. 
TC9 consists of a small core of dense gas and dust. The 
central source is surrounded by sufficiently large amount of 
materials that the 9.7/im silicate band is still detected in 
absorption. Besides, the SED of the central source is typi- 
cal of Class II objects, consistent with the active accretion 
phase being now terminated. We adopt as a rough estimate 
of the ionizing flux impinging on TC9 the one measured at 
the surface of TC2 : $ ~ lO^^cm-^s-i. This is justified 
because TC2 and TC9 have similar projected distances to 
the exciting star of the nebula. Using Eq. (36) of Lefloch & 
Lazareff (1994), which relates the mass-loss rate to the size 
of the condensation and the ionizing photon flux, we find 
a mass loss rate of ~ 8 Mq yr~^, where we have adopted a 
mean radius of ~ 8"; hence, the condensation will have been 
photoevaporated in a few 10^ yr, while the young stellar ob- 
jects TCI and TC2 will still be shielded from the ionizing 
radiation. TC9 is a clear case of a protostellar core overrun 
by the ionization front, which is about to disappear in the 
nebular gas. 

IRS 6-7, TCll-12,13 illustrate the ultimate stages of 
early stellar evolution when the parental envelope of the 
newly born stars has been almost fully photoevaporated 

and hardly any cold molecular gas emission is detected. In 
these objects, self-gravity fails to bind the structure, as ev- 



idenced by the ratio between the virial and the actual mass 
of the condensations. We note however that the gravita- 
tional potential due to the central object(s) was ignored, so 
that the actual discrepancy might be less. 

7.3. Core formation and Filaments 

The overall distribution of molecular material differs a lot 
from the free-free emission distribution observed with the 
VLA, which traces an almost perfect Stromgren sphere 
(CL98). Most of the protostellar cores are encountered in 
filamentary structures. Some of them, hke TCI, TC2, TC4, 
are detected in the large-scale photon-dominated region but 
its not the case for all of them. In particular, cores and frag- 
ments in WF, such as TCOO, 6, 7 are encountered far away 
from the HII region. 

Since the WF filament stretches over « 10 pc, far be- 
yond the ionization radius of M20 and connects several 
complexes together, we can exclude that it formed from the 
fragmentation of the layer of molecular material swept-up 
in the expansion of the HII region, as in the "collect and col- 
lapse" scenario (Elmegreen & Lada 1977; Elmegreen 1989). 
WF was already present before the ignition of the HII re- 
gion, and so were probably the dust lanes of the Trifid, and 
other filaments, like in the TC4 region. The separation of 
the massive condensations in WF along the filament is ap- 
proximately constant, equal to 1.1 x 10^^ cm~^ (3.6 pc); this 
provides additional evidence for the fragmentation scenario. 

In the filaments EL, SL, and TC4f, the size of the cores 
is very close to the diameter of these filaments, and their 
Jeans length. This fact supports the idea that the cores 
form as a consequence of the fragmentation of the filaments. 
Also, many of the pre-stellar cores, such as TC6,10,11 are 
not gravitationallybound, which implies that they probably 
formed dynamically inside the filaments. 

We explore here the scenario that all the condensations 
detected in our survey actually result from the fragmen- 
tation of filamentary clouds. These questions will be ad- 
dressed more thoroughly in a detailed analysis of the molec- 
ular gas dynamics in a forthcoming article (Lefloch 2008, 
in prep). 

7.4. Fragmentation of WF 

The map of the thermal dust emission shows that the mate- 
rial of WF is concentrated mainly in the massive condensa- 
tions TCOO, TCO, TC3, TC5. From the average gas velocity 
dispersion in the TC0-TC3-TC6 region {a = 0.7 km s'^), 
we estimate the virial mass per unit length m^ — 2(T^/G = 
220 A'/q pc~^, which compares well with the mass per unit 
length in WF (Sect. 4). We conclude that the filament was 
initially self-gravitating. Under the simplifying hypothesis 
that WF was initially in hydro dynamical equilibrium, and 
that it fragmented under some external perturbations, the 
properties of the fragments will be mainly determined by 
the fastest growing mode of the perturbation. 

We present here a simple analysis, based on the results 
of the continuum survey. The characteristic size for a fila- 
ment in hydrodynamical equilibrium is given by : 

H = 0.0035 ( ( — ^) pc 

V0.3kms-V \2x lO^cm-^/ ^ 
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where Cg is the effective sound speed (= Id velocity dis- 
persion a) and no is the gas density in the central region 
of the filament (see Nakamura, Hanawa & Nakano, 1993), 
and is approximately twice the average density of the fila- 
ment. Hence, we find H = O.lSpc and an effective filament 
radius of Re = 2\/2H = 0.51 pc. The actual diameter of 
the filament is not so easy to measure from the emission 
maps; the initial diameter of the filament could have been 
larger before undergoing external compression. The obser- 
vations at SEST and CSO show a somewhat larger radius 
(of about 0.6 — 0.8 pc), though not inconsistent with these 
results. For a purely hydrodynamically supported filament, 
the most unstable mode of axisymmetric perturbations has 
a wavelength A = 22.QH ~ 4.0 pc, a value in reasonable 
agreement with the mean separation between fragments 
« 400" = 3.6pc. 

A better agreement is obtained under the assumption 
that the filament is threaded by a magnetic field and that 
equipartition is reached between kinetic and magnetic en- 
ergies. The role of the magnetic field must be taken into 
account as it provides additional support against the fila- 
ment self-gravity; as a conscqiicncc the filament increases 
its scale height and its effective radius Ro. With the same 
physical parameters and assuming equipartition between 
magnetic and kinetic energy, the scale height Hi becomes 
H' = H X a/2 = 0.26 pc and the equivalent radius Re = 
0.73 pc (90"). The wavelength of the most mistable mode 
for axisymmetric perturbations in a magnetized filament is 
then A = 13.0if' = 3.4 pc. Overall, these results agree well 
with the observations, in particular the constraint on the 
diameter of the filament is better fulfilled. 

Inspection of the fragment along WF shows a marked 
trend. Whereas TC3 and TC5 fragments have very com- 
parable masses (« 800 M©), the mass of each newly en- 
countered fragment when moving North along the filament 
decreases, by a factor about 2 (see Table 2). The material 
distribution of TCOO seems less structured than in TOO, 
itself less structured than in TC3 and TC5. This suggests 
that the fragmentation of the filament could have started 
first with the formation of TC3 and TC5, then proceeded 
to the North of the filament, where TCO and TCOO would 
represent subsequent generations of fragments. 

The sequential star formation triggered by the expan- 
sion of an HII region, impacting a filamentary cloud has 
been studied numerically by Fukuda & Hanawa (2000). 
Their simulations show that under a wide range of physical 
conditions it takes about 10 times the radial sound crossing 
time 

T = H/Cs 

/ \ -1/2 

= 1.7xl0M 3) yr 

V2xl06cm-V ^ 

to form the first condensations in the shocked filament. In 
the case of WF, this time would be as large as 4 Myr, almost 
one order of magnitude larger than the age estimated for 
the Trifid (L02). Despite all the uncertainties in the age de- 
terminations, the detection of photoionized protoplanetary 
disks in the central cluster of the HII region (LOl) implies 
that the nebula is young. It is more likely that fragmen- 
tation in WF was initiated by some other events than the 
interaction with M20. An interesting result of these simula- 
tions is that fragmentation can occur at much earlier times 
as more kinetic energy is injected in the shock. An alternate 



Fig. 16. Map of the 1300/^m thermal dust emission super- 
posed on a map of the 20cm free-free emission of the region 
in greyscale (from Yusef-Zadeh et al. 2000). 



scenario is that fragmentation was already ongoing before 
the ignition of the nebula. 

A large number of observations support the existence 
of a physical interaction between the young W28 SNR 
and the surrounding molecular gas (Wootten 1981, Frail 
& Mitchell 1998, Arikawa et al. 1999; Rho & Borkowski 
2002; Yusef-Zadeh et al. 2000; Reach, Rho & Jarrett). It is 
not clear however if such events as energetic as supernova 
remnants could trigger the fragmentation of the filament on 
such short timescales. On the other hand, the characteris- 
tic timescales of TC3 and TC5 are fully compatible with a 
dynamical interaction between WF and W28. As discussed 
above, the free-fall gravitational timescale in the cores of 
TC3 and TC5 are extremely short, of the order 10^ yr, and 
similar to the age of the protostellar sources in TC3 and 
TC5. This is in excellent agreement with the estimated age 
of the W28 supernova remnant (3.3 x 10* yr, Velazquez et 
al. 2002). 

In order to confirm this hypothesis, it should be inves- 
tigated whether the impact of the SNR shock on the WF 
filament, the shock could trigger the gravitational collapse 
of the core and the formation of the protostellar cluster. 
The location of TC5 at the edge of the S28 shell star for- 
mation strongly suggests that star formation has indeed 
been triggered in the core as soon as the supernova shock 
hit the condensation. 



7.5. Fragmentation of TC4f 

In TC4f filament, the different cores (4a, 4b, 4c, 4d) are lo- 
cated along the filament, at the border of the nebula. A de- 
tailed analysis of the physical conditions showed evidence 
for compression of the filament gas (LCOO). No evidence of 
velocity gradient was found between the various molecular 
gas tracers. If the collapse of the cores had started before 
compression of the surface layers occurred, lines with dif- 
ferent optical depths should probe different regions of the 
core, hence tracing regions at different velocity. There is 
no such velocity gradient observed, between the optically 
thick CS and the optically thin C^^O line. Most likely, 
the cores TC4a-b-c-d formed after compression of the fila- 
ment. The properties of the cores are consistent with them 
being fragments of this layer (LCOO). The timescale for 
gravitational collapse of a decelerating shocked gas layer 
is T ~ 0.25(0/90)"^^^, where po is the typical (pre-shock) 
density of the parent molecular cloud (Elmegreen, 1989). 
The average density determined over WF (2000 cm~^) pro- 
vides a good estimate of po- The typical timescale for the 
fragmentation of TC4f density is therefore 0.35 Myr, com- 
patible with the age estimated for M20. To summarize, 
there is good observational evidence that TC4f is a fila- 
ment compressed by the expansion of the HII region, which 
subsequently became gravitationally unstable and formed 
young stellar objects. 
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7.6. Fragmentation in the dust lanes 

We have carried out a similar analysis for filaments SL and 
EL. The results are summarized in Table 3. Unlike WF and 
TC4f, the Southern and the Eastern lanes are only weakly 
gravitationally boimd. They arc probably confined by the 
external pressure of the ambient medium or by a helicoidal 
and/or toroidal magnetic field (see Fiege & Pudritz, 2000a). 
This is consistent with the less efficient star-forming effi- 
ciency observed there. However, the actual mass per unit 
length of the lanes is less since the O star turned on and has 
started to photoevaporate the surrounding molecular gas. 
It is not possible to accurately estimate the mass loss rate of 
both filaments as the orientation of the dust lanes makes it 
difficult to discriminate the contributions of the photoion- 
ized lane surface layers and the nebular gas to the free-free 
emission; moreover we do not know the exact distance of 
the O star to the lanes. A reasonable order of magnitude 
can be obtained by adopting the expression of the mass- 
loss rate derived for cometary globules to the filaments, 
and assuming an impinging ionizing photon flux similar to 
the one measured towards TC2. Under such hypothesis, we 
estimate a typical mass loss rate of 85 M© Myr~^ for a fila- 
ment 1 pc long. With an estimated age of 0.4 Myr, it means 
that SL has evaporated about half its initial mass, which 
was « 85 M© pc^^. Similarly EL has lost about 1/3 of its 
initial mass. In both cases, the initial mass of the filament 
was noticeably closer to the virial mass; hence, gravity was 
playing a more important role in the overall mechanical 
equilibrium of the clouds, which means that the discrep- 
ancy between filament and cores with respect to self gravity 
was probably much less. 

Fiege & Pudritz (2000a, b) have investigated the sta- 
bility and fragmentation of filamentary clouds truncated 
by an external pressure and threaded by a helical mag- 
netic field. They addressed the stability of the filaments to 
gravitational fragmentation and axisymmetric MHD driven 
instabilities. The authors have identified two well sepa- 
rated types of unstable modes, depending on the toroidal 
to poloidal magnetic fiux ratio F^/F^, which join together 
when r$/F^ = 2. At low values, the instability mode is 
gravity-driven and for pressure truncated filaments such as 
EL and SL, the growth timescale is ~ 1 Myr or longer. It 
implies that fragmentation would have started well before 
the birth of the exciting star of the nebula. The unstable 
MHD-driven instabilities triggered above F^/F^ = 2 can 
have very large growth rates compared with gravity-driven 
modes, and hence much shorter fragmentation timescales. 

The regime of instability relevant to the lanes can be 
determined from the ratio of the fragment separation A to 
the filament diameter D. In our case, one has X/D = 1.7 
for EL and 1.5 for SL, respectively. We note however that 
the photoevaporation of the outer, low-density, layers of the 
filaments underestimate the actual ratio, before ionization 
turns on. The radius of the filaments in the neutral parental 
cloud, was definitely larger, so that the ratio is probably 
closer to 1 or even less. Prom Pig. 17 of PPb it comes that 
F^/F^ = 2 — 3, i.e. MHD driven instabilities can propagate 
in the two filaments. The most unstable modes are charac- 
terized by a growth rate « (0.003 - 0.1) x (47rGpo)~^- 
The growth timescale of these modes (w~^) can be as low 
as 4 X 10^ yr in the case (F<i>/Fz ~ 3). This is of the same 
order as the radial soimd crossing time Re/cs = 0.3 Myr, 
a constraint to be met for an instability to develop along 



the filament without destroying the entire structure. We 
conclude that axisymmetric MHD-driven instabilities can 
in principle account for the fragmentation observed in the 
dust lanes of the nebula, on a timescale compatible with 
its age. However it puts some constraints on the magnetic 
field structure of the filaments. 

At this stage, our work suggests that the magnetic field 
probably plays an effective role in the dynamics of the fil- 
aments of the Trifid and in their fragmentation. Whatever 
the relative importance of hydrodynamical and magnetic 
support, the main conclusion is that the various condensa- 
tions TCOO, TCO, TC3, TC5 most likely formed from the 
fragmentation of WP. Both fast and slow fragmentation 
processes are at work in the Trifid, and depend on the ini- 
tial conditions in the environment of the nebula. Zeeman 
and polarimetric measurements of the magnetic field are 
needed to quantify its role in the support against gravity 
and the formation process (es) of the stellar incubators. 

8. Conclusions 

Our mapping of the Trifid Nebula at millimeter, submil- 
limeter and mid-infrared wavelengths has allowed us to 
characterize the gas and dust distribution. Thirty three 
cores have been identified, which cover all the stages of 
protostellar evolution, from the early prestellar to the late 
protoplanetary phases. Sixteen cores are found to display 
evidence of star formation. 

The cores are found to lie at various stages of photoion- 
ization, from the deeply embedded molecular phase to the 
cometary phase and the late stages of photoevaporation. 
Taking into account also the previously studied objects 
like TC2 (LC02), all cometary globules display very sim- 
ilar properties. They have typical sizes of 0.1 to 0.2 pc, H2 
densities ranging from « 10^ cm~^ to 10^ cm""^, and masses 
of a few Mq up to a few tens M©. They appear to form 
at most one intermediate-mass object. Hence, we conclude 
that photoionization does not increase the star formation 
efficiency in the individual cores. We find evidence of clus- 
tering only in the shocked, dense condensations, such as 
TC3 and TC4 (previously studied in LCOO) and TC5. 

Comparison of the line and continuum emission prop- 
erties shows hardly any evidence for molecular depletion 
inside the condensations. Only moderate molecular deple- 
tion has been reported towards the densest sources (TC3; 
LCOO). The mid-infrared 5 — 17/Ltm fiux detected with 
ISOCAM towards TCI and TC9 is consistent with the 
emission of a hot central source, with dimensions typical 
of protoplanetary disks, which is absorbed by the parental 
envelope of cold dust and gas. ISOCAM spectroscopy has 
revealed the presence of large amounts of CO2 ices in 
cometary globule TCI, thereby unveiling the potentially 
rich composition of icy grain mantles in protostellar envi- 
ronments in the Trifid. As the ice composition reflects the 
various processes that take place on grains and the chemi- 
cal evolution of the protostellar envelope, mid-infrared ob- 
servations should allow us to investigate how the chemical 
composition of the circumstellar envelopes varies along with 
protostellar evolution in the HII region. 

The cloud material is distributed in filamentary struc- 
tures of relatively high density (a few 10^ cm~'^) and gas 
column density N(H) of a few 10^^ cm~^. Their properties 
compare well with those of IRDCs and the Orion filaments 
reported by Maddalena et al. (1986). The length of these 
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structures ranges from about 1 pc to 15 pc or more. Since 
some of them extend far beyond the ionization front of the 
HII region, these filamentary structures cannot have formed 
in a "collect and collapse" manner. They were pre-existing 
the birth of M20 and reflect the initial conditions of the 
parent molecular cloud. 

Among the four filaments studied in greater detail by us, 
two of them are confined by self-gravity (WF and TC4f), 
whereas the other two are presently not. Taking into ac- 
count the photoevaporation of the surface layers of the lat- 
ter, they were probably self-gravitating too before the on- 
set of photoionization. All these filaments are fragmented 
in condensations with typical masses of a few hundreds of 
Mq . The fragmentation of the filaments can be accounted 
for by MHD-driven instabilities, although the exact impor- 
tance of the magnetic field in the region remains to be 
quantified. Comparison with numerical simulations shows 
that the impact of M20 on the filaments is too recent to 
have triggered their fragmentation, with the exception of 
the TC4f filament. It may well be, however, that this frag- 
mentation occurred at an earlier stage or was accelerated 
under the influence of high energy shock interactions with 
young supernova remnants. It is not clear if the mass distri- 
bution of the fragments in the filament WF results from the 
propagation of instabilities inside the filament, hence trig- 
gering the formation of the subsequent generations of cores. 
Conversely, the small-scale fragmentation reported in the 
filament TC4f provides a convincing example of a filament 
compressed and destabilized in the expansion of the HIT 
region, well accounted for by the modelling of Elmegreen 
(1989). The mass of the individual fragments formed in the 
TC4f filament is rather low (~ 10 Mq), whereas the mass 
of the objects in WF is typically ten times higher. 

The star formation going on in the dense cores TC3 and 
TC5 is characterized by very short dynamical timescales 
(about 10** yr), comparable to the age of the nearby W28 
supernova remnant. TC5 is a good candidate to investi- 
gate the possibility of SNR shocks as a trigger of star for- 
mation. Therefore, the environmental conditions appear to 
have played a very important role in determining the star 
forming conditions in the Trifid Nebula. The kinematics 
of the gas will be investigated in a forthcoming paper, in 
order to search for signatures of these environmental condi- 
tions and of possible interactions with the highly energetic 
phenomena which have recently taken place in the neigh- 
bourhood. 
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Appendix A: Cores in the Western Filament 

A.l. TC5 

TC5 is the brightest source discovered in our survey, with 
a peak flux of ~ lJy/11" beam. It is located in the south- 
ernmost fragment (« 810 Mq) of the Western filament WF 
(Figs. [TI IA.ip . The isocontours of the 20cm free-free emis- 
sion show that the core lies at the edge of the W28 super- 
nova remnant. 



Fig. A. 2. Molecular gas emission obtained towards TC5 in 
the millimeter transitions of SiO and CS. Fluxes are ex- 
pressed in units of antenna temperature. 



The core is elliptical with major and minor axis of 
0.20 pc X 0.13 pc, respectively. The brightness peak is off- 
centered, suggesting the presence of multiple sources in- 
side the continuum core. Indeed, SPITZER unveiled a clus- 
ter of sources in the dust central core (see right panel in 
Fig. IA.1|) . The emission is strongly concentrated towards 
the peak. A mass of 100 M© (about half the core mass) 
is contained within a 9" region, which indicates H2 densi- 
ties n(Il2) ~ 8 X lO^cm"^). The continuum source is not 
detected in the near- infrared with 2MASS. It however coin- 
cides with the IRAS pointsource 17589-2312. Comparison 
of the bolometric luminosity (4000 Lq) to the millimeter lu- 
minosity yields a ratio Lboi/(10"^ x L1.25) — 15. This value is 
low and typical of those encountered in Class protostellar 
cores. This testifies to the youth of the TC5 core. 

The line profiles display broad wings tracing gas that 
reach velocity of 20kms~^ with respect to the core, 
which testifies to the ongoing star formation at work in 
the core (Fig. IA.2[) . The discovery of methanol maser 
emission confirms indeed that massive protostars or non- 
ionizing intermediate-mass protostars are forming in TC5 
(Pestalozzi et al. 2005). Complementary interferomctric ob- 
servations are required to identify which of the continuum 
sources are powering the molecular outflows detected to- 
wards the core (Lefloch et al. 2008). The physical prop- 
erties of TC5 are therefore very similar to those of TC3 
and TC4. At the high densities encountered in the cen- 
tral core of TC5, the typical gravitational free-fall time is 
T = (37r/32Gp)i/2 ~ 1.2 x 10"* yr. The entire star-forming 
process taking place in the core is very recent. 

A.2. TC6 

The thermal dust emission shows a weakly contrasted 
condensation, without any dominant central condensation 
fFig. lA.3"|) . Instead, there are several components barely re- 
solved by the telescope beam, separated by typically 20". 
These sub-condensations (TC6A-D) have low masses of 
typically 2 - 3Mq. 

Only a few molecular lines were observed towards TC6, 
at SEST and CSO. As the telescope main beam has a 
size comparable to the dimensions of the condensations, 
the main-beam filling factor must be taken into account 
to accurately estimate the brightness of the CS lines and 
the density of the condensations. The filling factor / is 
0.81 and 0.65 for the J = 3 ^ 2 and the J = 2 ^ 1 tran- 
sitions, respectively, assuming a Gaussian source distri- 
bution / = O'^/{0'^ + 9%). Hence, the CS line brightness 
temperature of the J = 3 ^ 2 line components are 0.49 K 
and 0.35 K, respectively, and 0.74 K and 0.65 K for the 
J 2 ^ 1 line. 

Neither ISO nor SPITZER detected pointhke sources 
towards TC6, which is consistent with TC6 being in a 
prestellar phase (see Fig. [T|). Indeed, the mass of the sub- 
condensations is much less than the virial equilibrium mass 
(see Sect. 7). Two kinematic components are detected in 
the CS lines (see Table 5). The line profiles are broader 
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Fig. A.l. Contour map of thermal dust emission at 1.3mm observed with the IRAM 30m telescope superposed on an 
Sum image of the region (greyscale). Contour levels are 25, 50, 75, 100, 150, 300, 450, 600, 750,900 mJy/11" beam. The 
Western filament WF is detected in absorption against the background. First contour and contour interval are 5 mJy/11" 
beam. A magnified view of the protostellar cluster associated with the TC5 core is presented in the right panel. 



Fig. A. 3. Thermal dust emission map of TC6 at 1.3mm 
observed with the IRAM 30m telescope. First contour and 
contour interval are 5mJy/ll" beam. 



Fig. B.l. Molecular line spectra obtained at the emission 
peak of the molecular core TC13. Fluxes are expressed in 
units of antenna temperature. The dashed line marks the 
velocity of the main body gas. 



Fig. A. 4. Millimeter line emission observed at SEST in 
TC6 at the offset position (-280", -440"). The CS Unc in- 
tensities have been corrected for main beam dilution effects. 



than those observed in the dust lanes of the nebula, and 
indicate a higher degree of turbulence. 



Appendix B: Cores in the Eastern Lane 

B.l. TC8 

Figure [7] shows a second condensation gas in the Eastern 
lane. Emission of the dense layers gas was observed in CS 
in the J = 2 ^ 1 transition at IRAM and in the J = 3 2 
transition at SEST. Unfortunately, the fragment and TC8 
could not be fully mapped. 

The molecular emission peaks near the offset position 
(+93. 6", +88. 6") where a Class I protostar has been dis- 
covered by SPITZER (Rho et al. 2006). The Unes of CO 
and its isotopologucs have intensities comparable to those 
measured in TCI (Fig.[Tni), indicating a kinetic temperature 
of about 30 K. Due to incomplete mapping, we have esti- 
mated the physical conditions in the core from the emission 
at the position (+90", +80") close to the SPITZER contin- 
uum source. Because the core size is of the same order as 
the telescope beam size, the effects of the convolution with 
the main beam cannot be neglected, unlike the case of TCI, 
whose extent is typically twice the telescope beam size. The 
intensity of the CO J = 2 — 1 line is « 0.75 K at the bright- 
ness peak. The beam filling factor reads 0^/(9^ +9%) where 
9s = 20" and 9b — 25.3" are the source and beam size 
(HPFW), respectively. A good estimate of the line bright- 
ness is obtained by dividing the main-beam brightness tem- 
perature by the beam filling factor. In the case of the CS 
J = 2 1 line, the beam filling factor is ~ 0.385, and the 
brightness temperature of the line is 1.9 K. From the main- 
beam brightness temperature measured at SEST, we esti- 
mate a brightness temperature of 2.2 K for the J = 3 ^ 2 
line (the beam size is 33.7" at this frequency). Hence, the 
line intensities compare well with those measured towards 
TCI. 

As a conclusion, TC8 is characterized by molecular gas 
and dust emission which is less constrated than in TCI. 
However, the mass and the physical conditions in both con- 
densations are rather similar. The dense gas region in TC8, 
as traced by CS, occupies a smaller part of the molecular 
gas fragment. The mass of the CS core represents only 25% 
of the fragment mass associated with TC8. In TCI, the CS 
core represents 80% of the whole mass. We interpret this 



Fig. B.2. Molecular gas emission observed towards TC13, 
integrated between 2 and 8 km s^^. The arrow points to the 
position of the protostellar source. 



Fig. C.l. Molecular line emission in the Southern Lane re- 
gion, integrated between -3 and +3kms~^, showing TCIO 
and TCll. 



Fig. C.2. Molecular line spectra observed at the emission 
peak of TCIO in the Southern Lane. Fluxes are expressed 
in units of antenna temperature. 



difference as a result of a more effective photoionization and 
compression of TCI. 

B.2. TC13 

The Eastern lane bends in the Southeastern direction be- 
yond TC8 (see Fig. [T]). This part of the dark lane contains 
lower amounts of gas and dust. As a consequence, it is de- 
tected in molecular gas emission (CO isotopomers and CS) 
while it is missed by our 1.3 mm continuum mapping. A 
condensation of 55" x 25" was discovered in this part of the 
filament (Fig. IBT2| . A Class O/I source candidate (TC13) 
was identified by SPITZER, which lies in the outer layers 
of the condensation. 

The CS J 2 ^ 1 transition was detected at SEST. 
The line has a typical main-beam brightness intensity 
Tmb = 0.6 K. Once corrected for beam dilution (the beam 
size is comparable to the size of the core), we estimate 
a brightness temperature — 1.3 K at the brightness 
peak (Fig. IB.1[) . Assuming a standard relative abundance 
[CS]/[H2] = 10-^ the CS column density is N(CS) ~ 
6 X 10^^ cm^^. The gas density must be ~ 10^ cm~'^ in 
order to account for the brightness of the CS J = 2 ^ 1 
line. 



Appendix C: Cores in the Southern Region 

Two condensations, TCIO and TCll respectively, could be 
identified in our molecular emission maps (see Fig. IC.ip . 
The two condensations are separated by « 84" = 0.7 pc. 
Only the first condensation is detected in the millimeter 
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continuum. The overall emission is weak; the peak flux is 
about 25 niJy/11" beam, and the integrated flux is 25 niJy, 
implying a peak column density N(H) = 2.3 x 10^^ cm~^. 
The central core, defined by the contour at half power, has a 
typical size of 7" (beam-deconvolved) and a mass of 2.5 Mq. 

The molecular gas condensation, defined by the contour 
at half-power, has an elliptical shape with major (minor) 
axis of 80" (30") contains a mass of 25 M0, almost one 
third of the material contained in the lane (66 Mq). 

Lines are somewhat narrower in the Southern lane. 
We measure T^i^ = 2.5 K (Au = 1.4 km s^^) (Fig. [C2I. 
Assuming a standard relative abundance [CS]/[H2] = 10~^, 
the CS line emission is consistent with H2 densities of 
2.5 X 10^ cm^'^. In the Northern region, the molecular gas 
emission reveals TCll as a larger condensation of larger di- 
mensions and lower gas column density. Adopting the same 
physical conditions as in the South, we derive a peak col- 
umn density N(H2) = 9 x 10^^ cm~^ and a mass of 6 Mq 
from the ^^CO data. No evidence of protostellar activity 
was detected in the Southern lane. This is consistent with 
the mass of TCIO and TCll being much lower than their 
virial mass. 
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